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Abstract 



We have compared the virtual corrections to single hard bremsstrahlung as 
calculated by S. Jadach, M. Melles, B.F.L. Ward and S.A. Yost to several 
other expressions. The most recent of these comparisons is to the leptonic 
tensor calculated by J.H. Kiihn and G. Rodrigo for radiative return. Agree- 
ment is found to within 10~ 5 or better, as a fraction of the Born cross section, 
for most of the range of photon energies. The massless limits have been shown 
to agree analytically to NLL order. 



f Work partly supported by the US Department of Energy Contract DE-FG05- 
91ER40627 and by NATO grant PST.CLG.980342. 



High precision studies of the Standard Model at proposed linear colliders 
will require per-mil level control of both the theoretical and experimental 
uncertainties in many critical processes to be measured. This will require 
computing higher order electroweak radiative corrections at least to order 
0(a 3 L 3 ) for the leading log effects and to the exact 0(a 2 ). One important 
contributions at this order is the virtual photon correction to the single hard 
bremsstrahlung in e + e~ annihilations [1-3,6,7]. 



Figure 1: Feynman graphs for the virtual 0(a) correction to the process 
e+e" -» fj + 7. 

In Refs. [3], some of us (S.J., B.F.L.W., S.A.Y.) have presented compar- 
isons of the results in Refs. [1-3], and in general a very good agreement was 
found. In particular, all of these results can be shown to agree analytically 
to NLL order. However, some differences become apparent at the level of the 
NNLL (next-to-next-to leading log), depending on the different levels of "ex- 
actness" in the calculations. Specifically, the mass corrections are included in 
Ref. [3] in a fully differential way, whereas in Ref. [1] the mass corrections are 
included but the photon angular variable is integrated over and in Ref. [2] the 
results are fully differential but the mass corrections are incomplete. These 
comparisons therefore can not really test the NNLL, fully differential results 
in Ref. [3]. 

More recently, another fully differential expression for the virtual photon 
correction to single hard bremsstrahlung has appeared in Refs. [6,7], where 
the result is obtained by an independent method via a leptonic tensor cal- 
culated to investigate radiative return, with particular emphasis on the 1-2 
GeV cms energy regime. This permits a cross-check at the NNLL level of the 
corresponding results. Details have been presented in Refs. [4,5]. An indirect 
comparison of the results in Refs. [3, 6, 7] was also reported in Ref. [8] via 
a comparison of the two Monte Carlo programs PHOKHARA [9] and /C/C 
MC [10], as these two Monte Carlos have the realizations of the results in 
Ref. [6,7] (PHOKHARA) and Ref. [3] ( /C/C MC). Agreement at the per-mil 
level was found on selected observables. 
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Relevant Feynman graphs for the process e + e~ — > ff + 7 are illustrated 
in Fig. ^ in Ref. [3] the ISR matrix element is evaluated using the helicity 
spinor amplitude method of CALKUL [11], Xu et al. [12] and Kleiss and 
Stirling [13]. Computer algebra techniques [14] were used for reducing the 
loop integrals. 

The 0(a 2 ) virtual correction to single hard bremsstrahlung can be ex- 
pressed in terms of a form factor multiplying the 0(a) tree level matrix 
element [3]: 

^i SR(1) = £(/„ + fih + hh)M?^ (1) 

where A / li SR( ' ^ is the tree- level hard bremsstrahlung matrix element and 
.M* includes an additional virtual photon. The form factors f and spinor 
factors Ii can be found in Ref. [3] or Ref. [5]. (The latter form is equivalent 
to the earlier form, but is explicitly numerically stable in collinear limits.) In 

— (n\ 

the Monte Carlo program, we actually calculate the YFS residuals [16] (3\ 
which are obtained by subtracting the soft limit 47t5yfs defined in Ref. [3] 
from /o in Eq. ^ and using the matrix element to compute the cross-section. 

The mass corrections are then added following the methods in Ref. [15], 
after checking that the exact expression for the mass corrections differs from 
the result obtained by the latter methods by terms which vanish as m 2 e / s — > 
where m e is the electron mass and s is the squared cms energy. The effect is 
to add a correction 

it v 2m l f r i , r 2 

\J0/m — 



r 2 Ti) (1 - ri) 2 + (1 - r 2 ) 2 
x I (Jo) + (In z - l)L - Z - In z + ± In 2 z + l} (2) 

to the spin-averaged form factor (/)o, where L = ln(s/m 2 ) is the "big log- 
arithm" in the LL expansion, z = s'/s and = 2pi ■ k/s, with incoming 
fermion momenta Pi and real photon momentum k defined as in Ref. [3] . 

To NLL order, an expression for the cross section which correctly repro- 
duces all collinear limits can be obtained by setting fi — f 2 = and using 
the spin-averaged NLL form factor 

+ 21nr 2 ln(l - r x ) - ln 2 (l - n) - ln 2 (l - r 2 ) + 31n(l - n) 

+ 31n(l-r 2 ) + 2Sp(r 1 )+2Sp(r 2 ) + (/ )r L (3) 
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where Sp(x) is the Spence dilogarithm function and only terms surviving 
when 7"i — > or r 2 — > are needed in the mass correction. In the NLL limit, 
the 0(a 2 ) YFS residual (3 1 for real + virtual emission may be expressed in 
terms of the YFS residual j3 x X for pure real emission using the result Eq. |HJ 

Jl ,, ^,"(i^(/.r). (4) 

The new comparison uses the leptonic tensor of Ref . [7] , 



Lq V = —pi {avaS-qT + anPiPi + (I22P2P2 + a 22 + P1P2) 
ss' 

+ i7ra_i(pfp2+Pi3^)} ( 5 ) 



The coefficients may be found in Ref. [7] The squared matrix element for 
real photon emission can be obtained by contracting this with a final-state 
tensor 

= e 2 (pftfl + vlvl - p 3 ■ ptrT) (6) 

to obtain 

|.M ISR | 2 = zL^H, u . (7) 

The coefficients can be decomposed as = a$ + ^ a\f in terms of 
a real-photon emission term and a virtual correction of order a. The YFS 

residual (3 1 may be obtained by calculating the cross section using Eq. [7| 
with a$ replaced by = — a* R , where the IR-divergent term is 



a m = a (0) 
13 »j 



1 7T 2 

2(L-l)lnw min + -L-l + y 



(8) 



with cutoff v min on the fraction of the beam energy radiated into the real 
photon, v = 2E 1 /y/s. Again, it is possible to find a relatively compact 
expression for the NLL contribution to the YFS residual in the massless 
limit. The result is 

M 2) ^fi + ^cO+q, 0) 



2tt 



with coefficient functions 



C\ = ~7ri\ ( — + ZC22 - 2Ci2 



2a 00 v z 

Cn £C 2 2 C12 , , 

C 2 = ^ + — - T -coo, (10) 



where the IR-finite coefficients are to be evaluated in the collinear limits 
r\ — > or r 2 — > 0. We have verified that in these limits, C\ = (/o) NLL an d 
C*2 = 0, so that the results agree analytically in the massless NLL limit. 
The results are illustrated in Fig. |21 for the case // = We show 

— (2) 

the complete f3\ distribution for our exact result JMWY as presented in 
Ref. [3], both with and without mass corrections, the result IN of Igarashi 
and Nakazawa et al. [2], the result BVNB of Berends et al. [1], and the new 
comparison with exact results KR of Kuhn and Rodrigo in Ref. [7] (with 
mass corrections) and Ref. [6] (without mass corrections). The results were 
obtained using the YFS3ff MC generator (the EEX3 option of the K.K. MC 
in Ref. [10]) with 10 8 events. The NLL contribution calculated in Ref. [3] 
has been subtracted in each case to permit a clear NNLL comparison. We 
see that there is a very good general agreement between all of these results. 

Specifically, of the results agree to within 0.4 x 10 -5 for cuts below 0.75. 
For cuts between 0.75 and .95, the results agree to within 0.5 x 10~ 5 , if 
the result of Ref. [1] is not included, and within 1.1 x 10~ 5 if that result 
is included. For the last data point, at v — 0.975, the result of Ref. [1] 
is approximately 1 x 10 -4 greater than the others and is off-scale in Fig. |21 
while the remaining results agree to 3 x 10~ 5 . The difference between the KR 
and JMWY result attributable to differences in the mass corrections never 
exceeds 10 -5 over the entire range of u max - 

These results are consistent with a total precision tag of 1.5 x 10~ 5 for 
our 0(a 2 ) correction f3[ 2 ^ for an energy cut below v = 0.95. The NLL effect 
alone is adequate to within 1.5 x 10~ 5 for cuts below 0.95. The NLL effect 
has already been implemented in the /C/C MC in Ref. [10] and the attendant 
version of KK MC will be available in the near future [17]. 

These comparisons show that we now have a firm handle on the preci- 
sion tag for an important part of the complete 0(a 2 ) corrections to the // 
production process needed for precision studies of such processes in the final 
LEP2 data analysis, in the radiative return at $ and B-Factories, and in the 
future TESLA/ILC physics. 
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Figure 2: NNLL contribution J3{ — /9{ NLL for 10 8 events generated by the 
YFS3ff Monte Carlo as a function of the cut w max on the fraction of the beam 
energy radiated to the photon. The results is in units of the non-radiative 
Born cross section for e + e~ — > fi + fi~. 
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